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DEFINITION OF THE DESIGN TRAJECTORY AND ENTRY FLIGHT 
CORRIDOR FOR THE NASA ORION EXPLORATION MISSION 1 
ENTRY TRAJECTORY USING AN INTEGRATED APPROACH AND 

OPTIMIZATION 


Luke W. McNamara* * Robert D. Braun] 


One of the key design objectives of NASA’s Orion Exploration Mission 1 (EM- 
1) is to execute a guided entry trajectory demonstrating GN&C capability. The 
focus of this paper is defining the flyable entry corridor for EM-1 taking into ac- 
count multiple subsystem constraints such as complex aerothermal heating con- 
straints, aerothermal heating objectives, landing accuracy constraints, structural 
load limits, Human-System-Integration-Requirements, Service Module debris dis- 
posal limits and other flight test objectives. During the EM-1 Design Analysis Cy- 
cle 1 design challenges came up that made defining the flyable entry corridor for 
the EM-1 mission critical to mission success. This document details the optimiza- 
tion techniques that were explored to use with the 6-DOF ANTARES simulation 
to assist in defining the design entry interface state and entry corridor with respect 
to key flight test constraints and objectives. 


INTRODUCTION 

The NASA Orion Exploration Mission 1 (EM-1) is scheduled to be launched as early as 2017 
from the Kennedy Space Center. This is scheduled to be the next Orion flight following the com- 
pletion of Orion Exploration Flight Test l 2 . This uncrewed EM-1 mission will provide a critical 
step forward in preparing Orion for future missions. The Orion spacecraft, assembled with a Ser- 
vice Module provided by the European Space Agency, is scheduled to be launched using the newly 
designed NASA Space Launch System (SLS). After a successful EM-1, EM-2 will become the first 
crewed mission to use the Orion spacecraft and SLS. 

The EM- 1 mission will be a mission 25 day mission to a Distant Retrograde Orbit (DRO) and will 
spend approximately 5 days in the DRO. After the DRO the spacecraft will return to Earth and re- 
enter the atmosphere at approximately 11 ^ (36,090 ^). Once inside the atmosphere, the vehicle 
will fly a guided entry using Guidance, Navigation, and Control (GN&C) bank angle modulation. 
The vehicle will target a splashdown location in the Pacific Ocean west of San Diego, CA. 

The GN&C system works to enable the completion of several mission objectives. The navigation 
subsystem provides the guidance subsystem with deduced acceleration, velocity, and position in 
translation and rotation. The guidance subsystem uses that information to target the desired landing 
site using the guidance algorithm PredGuid 8,9,10,11 . PredGuid pairs a modern Numeric Predictor 
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Corrector (NPC) with the flight-tested Apollo Final Phase guidance algorithm. 12 This combination 
is designed to perform the skip-entry at Earth which enables the vehicle to reach the target landing 
site from greater ranges. To reach the target landing location the guidance algorithm calculates the 
necessary bank angle to decrease the expected error in touchdown range-to-target throughout the 
flight. 

Entry Sequence 

The entry sequence will begin at the separation of the Crew Module (CM) from the Service 
Module (SM) prior to the Entry Interface (El). The geodetic altitude of El where the atmospheric 
density of Earth begins to have an impact on the trajectory is 121,920 m (400,000 ft). Between 
CM-SM separation and El the vehicle is rotated heat shield forward to prepare for entry. During 
the hypersonic, supersonic, and subsonic portion of the trajectory before the Landing Recovery 
(LRS) sequence is initiated, the vehicle is banked by guidance pointing the vehicle lift vector in 
order to target the landing site. During the LRS sequence, as currently planned, the Crew Module 
will deploy Forward Bay Cover (FBC) pilot parachutes pulling away the FBC, deploy two drogue 
parachutes, release the drogue parachutes, and deploy three pilot parachutes which pull out the three 
main parachutes. The described parachute sequence for the EM- 1 Entry, Descent & Landing (EDL) 
flight is illustrated in Figure 1 below. 
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Figure 1. Orion EMI: Nominal Parachute Sequence 
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Trajectory Generation 

Trajectoies were generated with Advanced NASA Technology Architecture for Exploration Stud- 
ies (ANTARES) 7 , a six degree-of-freedom simulation, that numerically propagates the state vectors 
of the simulation bodies forward in time. ANTARES was configured to use prototype flight soft- 
ware rather than Rapid Algorithm MATLAB/Simulink® Engineering Simulation (RAMSES) flight 
software to increase flexibility and decrease required Computer Processing Unit (CPU) run time. 
The PredGuid guidance algorithm that was employed uses a numerical predictor-corrector com- 
bined with the Apollo Final Phase Entry Guidance. The PredGuid guidance gain set and reference 
trajectories were retuned or at least recalculated for varying and entry range to keep the expected 
trajectory performance including uncertainties inside the energy bucket described by Bairstow. 10 
To analyze the impacts of uncertainties in the design process, Monte Carlo (MC) analyses are used 
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to model the uncertainties as random dispersions on the simulation models and inputs. For example 
Monte Carlo dispersions are applied to the aerodynamic characteristics, initial position vector and 
velocity vector, initial body attitude and rates, atmospheric characteristics, monthly atmosphere, 
initial navigation errors, Global Position System (GPS) errors, Inertial Measurement Unit (IMU) 
eiTors, aerodynamic Reaction Control System (RCS) jet interactions, RCS jet thrusts and specific 
impulses etc. 

The El position and velocity are dispersed using a 6x6 state covariance in the UVW frame also 
called the RSW frame. 16 The expected state covariance is progressively being refined with more 
analysis and testing. The atmospheric dispersions are dispersed using the Global Reference Atmo- 
sphere Model (GRAM) 2010 version 3.0. 13 

ANTARES was configured to simulate the trajectories using the physics models given in Table 1 
and GN&C models in Table 2 among others. 


Table 1. Physics Models 


Physics 

Computer Model 

Atmosphere 

GRAM 2010 V3.0 13 

Aerodynamics 

MPCV Aerodynamic Database V0. 7 1.1 14,3 

Aero thermal Dynamics 

MPCV Aerothermal Database V0. 93 

Gravitation 

NASA’s JEOD Version 1.5.2 


Table 2. GN&C Models 

GN&C Component 

Computer Model 

Guidance 

Navigation 

Control 

PredGuid C-Code Flight Software 

Discreet Extended Kalman Filter using GPS, Barometric Altimeter, IMUs 
Proportional-Derivative Control using RCS 


ENTRY FLIGHT CORRIDOR DESIGN 

The capability of an entry vehicle to accomplish its objectives and meet its multitude of con- 
straints is a challenge. Typically constraints from other subsystems bound available domain space 
in most or all directions. 

Problem Description 

The purpose of defining the flight corridor is to capture and integrate constraints to understand 
the flyable domain of the vehicle. Once the constraints are integrated the corridor needs to protect 
for dispersions and facilitate design sensitivities 1 . 4 The current objective is to minimize CM mass 
at touchdown. To meet this many options are possible but the hope is to not increase operational 
complexity to accomplish it. As such, the GN&C objective is to find a suitable minimum value 
of lift-over-drag, j^*, to reduce the amount of ballast mass needed to position the center-of-gravity 
(CG) 6 accordingly. To achieve an increased jj value requires an increase in the required ballast 
mass. The is adequate if the design satisfies all the probabilistic constraints for predicted uncer- 
tainties modeled using Monte Carlo simulations. Such uncertainties are on atmospheric properties, 

‘For the meanings of variables refer to the Notation section near the end of the document. 
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initial position, initial velocity, reaction control system effective thrust, aerodynamics etc. A sec- 
ondary goal is to gain insight into the sensitivity impacts of independent variables and dependent 
variables on the integrated corridor. 

Hence, find the solution jj mrrL which minimizes the performance index F(.x) in Equation 1 where 
x is given in Equation 2, the available corridor width, c(x), is given in Equation 3 and the require- 
ment for protected corridor width, z, currently equals 0.15 degrees (previously was 0.30 degrees). 


F(x) = c — z 


( 1 ) 
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( 3 ) 


As the combined aerothermal and TPS structures groups desire to save TPS mass, and findings 
have indicated the result that this occurs at shallower flight-path-angles the objective from this group 
is 7 max of an available corridor. As such, that secondary objective is treated after finding the ^ 
from the GN&C design corridor. 

To define the design cycle for an entry corridor is subjective given the multiple independent 
variables. The technique employed here uses the following path of steps but there are many other 
paths that would lead to the same conclusion. 

First, define the mission goal to be to return from a specific orbital body or point in space defined 
for EM-1. That defines the El state return velocity range. Then use the minimum acceptable down- 
range to satisfy the crew HSIR, disposal concerns etc. Typically the minimum downrange case is 
defined to provide divert capability to avoid a large disruptive weather system such as a hurricane 
in the Pacific Ocean. Then follow the design cycle illustrated in Figure 2 below to converge on an 
acceptable corridor. 



Figure 2. One Version of a Design Cycle for Entry Corridor 
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If after the last step in Figure 2 the solution has been identified (converged) the process is com- 
plete. If not, then pick new heat load/rate limits and repeat steps 2 through 6. Then repeat the 
process for the longest beneficial downrange case. Due to the geometry of the Earth this range is 
8890 km (4800 n.mi) to have the ability to return to the desired landing site (based on its latitude) 
at any time during the lunar month. 

The key design constraints of the EM-1 entry trajectory corridor design are listed below with their 
probability success criteria. 

1 . Guided Vehicle Landing Accuracy. The landing accuracy requirement is for the vehicle to 
land within 10 km (5.4 nautical miles [n.mi.]) of the desired landing site with a probability of 
success of 99.865% (3er) at 90% confidence. As this requirement is treated in two pieces, a 
steep piece and a shallow piece, this is treated as two separate single tail distributions. Thus, 
this probability equates to allowing 1 miss for a 3000 case Monte Carlo analysis. 

2. Service Module Debris. The SM debris requirement is for potential SM debris to land out- 
side of buffer zones around foreign and domestic landmass with a probability of success of 
99.865% (3 <t) at 95% confidence. This probability equates roughly to allowing 0 failures 
for a 3000 case Monte Carlo analysis. This requirement is based on NASA Technical Stan- 
dard (NASA-STD 87 19.14), an accompaniment to the NPR for Limiting Orbital Debris (NPR 
8715. 6A), and is designed to limit the risk of human casualty due to end of mission disposal 
of jettisoned items. 

3. Guidance Degradation. It is desired that the entry guidance not succumb to atmospheric 
gradients that yield degraded performance. This protection ensures that guidance is able to 
protect against cases with HSIR violations between El and drogue parachute deployment. 

4. Guided Human System Integration Requirements (HSIR). The subset of Human System 
Integration Requirements (HSIR) this is concerned with are a set of requirements on the sus- 
tained crew member sensed linear acceleration, rotational acceleration, coupled rotational 
acceleration, rotational velocity and linear jerk. The requirement is to be satisfied for a prob- 
ability of success of 99.865% (3cr) at 99% confidence which equates roughly to 0 failures in 
a 3409 case sample size Monte Carlo. 

5. Ballistic Human System Integration Requirements (HSIR). The subset of Human System 
Integration Requirements (HSIR) this is concerned with are a set of requirements on the 
sustained crew member sensed linear acceleration, rotational acceleration, coupled rotational 
acceleration, rotational velocity and linear jerk. For the back-up mode of ballistic entry, 
executed in case of certain failures, the sensed linear acceleration limits are relaxed from the 
nominal guided entry case. The requirement is to be satisfied for a probability of success of 
99.865% (3<r) at 99% confidence which equates roughly to 0 failures in a 3409 case sample 
size Monte Carlo. 

6. Ballistic Entry Altitude Rate. The geodetic altitude rate is required to be less than or equal to 
zero to ensure successful atmospheric capture with a probability of success of 97.725% (2 o) 
at 90% confidence which equates roughly to 57 failures in a 3000 case sample size Monte 
Carlo. 

7. Ballistic Entry Structural Acceleration Load. The Crew Module’s structural acceleration 
load is required to be less than or equal to 15 g’s. The requirement is to be satisfied for 
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a probability of success of 97.725% (2<r) at 90% confidence, which equates roughly to 57 
failures in a 3000 case sample size Monte Carlo. 

8 . Aerothermal Thermal Protection System (TPS) Backshell Panel Maximum Tempera- 
ture. This aerothermal TPS requirement is currently booked at constraining the trajectory to 
a value less than or equal to 3150 °F. 

Listing all of these constraints in one equation yields 


' 91 ' 


Guided Landing Accuracy <10 km (Steep Side) 

92 


Guided Landing Accuracy < 10 km (Shallow Side) 

93 


Service Module Debris Disposal Limit 

94 


Ballistic HSIR 

95 


Ballistic Geodetic Altitude Rate< 0 

96 


Guidance Degradation 

97 


Guided HSIR 

98 


Ballistic Structure Acceleration < 15 g’s 

. 99 . 


Aerothermal Backtile Temperature < 3150 °F 


Problem Solution Subspaces. After multiple corridor design cycles, due to system level con- 
straints, the subspace of permissible independent variable combinations are defined as shown in 
Table 3 below. 


Table 3. Down Selected Design Space as a Result of Mass Reduction Activities 


Independent Variable 

Minimum 

Maximum 

j) (nd) 

0.23 

0.27 

El velocity, vj, (™) 

10,950 

11,050 

CM Mass, to, (kg) 

9,934 (21,900 lbm) 

10,387 (22,900 lbm) 

Downrange, R, (km) 

2,408 (1,300 n.mi.) 

6,482 (3,500 n.mi.) 

Flight Path Angle Corridor Width, c, (degree) 

0.15 

0.15 


DETERMINING CONSTRAINT CURVES 

In order to assess the safe flyable corridor of the subspace the constraints need to be defined 
numerically using ANTARES simulations. As knowledge of the constraint curves are critical in 
assessing the flyable corridor the curves would ideally be identified to a high degree of precision. 
This interest in increasing the curve fit of constraint lines needs to be balanced against the associated 
computational demands and wait time. The computational demands are not limited to CPU run time 
but rather also include disk memory to store the necessary data. 

The dependent variable constraint curve is defined by the associated independent variable set 
defined in Equation 2. Equation 2 can be rewritten into Equation 5 where range, R, has replaced 
longitude, A, and <5 has been removed (by tying the database to a fixed landing site). 


independent 


v„ 


1e 




R c 


171 El D Mach=25 


(5) 


Hence for a given Xindependent the dependent variable constraint curve could be characterized in 
all directions of the design space. To do that requires discretizing all dimensions of the Xindependent , 
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simulating a trajectory from each element of the discretized set in time from El to touchdown, and 
processing all of the dependent variable data to identify each constraint curve. To plot the constraint 
curve in a figure showing more dimensions than two is problematic for ease of understanding. As 
such the curve can be shown in any two dimensional domain such as jj Mach= 2 b vs - 7 b/ > Vm vs - 
7 ei, or jj Ma ,ch - 25 vs ' ^ei- Tis the intent is to locate the jj n/in that provides a r y EI corridor width of 
0.15 degrees it makes the most sense for this case to use the M ac .h= 2 5 vs - 7 ei domain to visualize 
the corridor. 

Hence for a certain set of mission re-entry constants 

T 

x constant = \ V E I 'f ei 77 EI m El ] ( 6 ) 

the entry comdor can be defined and visualized with respect to the independent variables 


Xindependent 


Iei 


L 

D Mach = 25 


1 T 


(7) 


and all the dependent variables can be found that define the constraints on the entry corridor in 
that two dimensional space. 

There are only two parts to constraint curve fitting to be concerned with here, finding (7, jj) 
points on the curve, and deciding the order of the polynomial to fit to the curve. Ideally those points 
on the curve would be from individual dispersed Monte Carlo sets and there would be a sufficient 
amount to characterize the curve in all directions of the design space. Currently that is unrealistic 
given the computation capabilities of the team and also it would be unnecessary considering the 
nature of the inputs to the problem from multiple subsystems. It is realistic however to characterize 
a curve in two dimensions, 7 and , either using undispersed trajectories for speed and breadth or 
using dispersed Monte Carlo trajectory sets for statistical accuracy. Several methods are described 
below that use variations of these methods. In addition, to find the points on the curve there are 
two sets of techniques, namely grid search using brute force predetermined constant increment 
parametric sweeps or iteratively using an optimization scheme. Although due to the discontinuities 
of the dependent variables and trip points, finding points iteratively using an optimization scheme 
presents certain pitfalls. 

Database of Undispersed Trajectories. A database of undispersed trajectories was generated for 
a discretized space of the same order as the number of independent variables in Equation 5, namely 
7. 


Xindependent, discretized 


y E l^\{Vmin Vmax\ 
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— D — D 7 


J J 
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Here the Xi t min and Xi, rn ax are defined wide enough to encapsulate any potential designs using 
input from orbit phase trajectory designers and the mass properties team. As such, for each element 
of the Xindependent, discretized set a multitude of dependent variables is available to extract constraint 
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lines. This database is used to create undispersed corridors and perform sensitivities to all the 
independent variables quickly. 

From Database of Undispersed Trajectories to Requirement Verification of Dispersed Trajec- 
tories. 

Now as the design requirements are based on certain probability of success and confidence val- 
ues including uncertainties, which are accounted for in trajectory design using dispersions on the 
trajectory models, an entry corridor would be more suitable by including the dependent variables 
(results) of the dispersed trajectories. 

To do this the corridor lines need to be redrawn through the targeted x constant that resulted in 
just meeting the requirement. The way this is done is using Monte Carlo. From binomial statistics 
the required probability of success and confidence level (value) are used to determine the necessary 
Monte Carlo sample size and number of cases from that sample size that are allowed to fail (not 
pass) the success metric. 

Using Monte Carlo Anchor Points 

Monte Carlo points identifying the constraint in the (7, |j) space are necessary to anchor the 
points with statistical meaning. Those points are found using an iterative method or with a grid 
searching method. In either case the constraint line is identified by analyzing adjacent sets of Monte 
Carlo’s and finding the tripping point where the dependent variable begins to surpass the constraint. 
In the iterative method an automatic setup using a simple optimizer such as a line search can run 
a Monte Carlo, assess all of the dependent variables to determine if the augmented root (constraint 



means iterating with values of 7 for a set constant , for every j-> anchor point needed, and repeating 

this process for every constraint. The convergence limit on the 7 iterations can be set according to a 
fidelity 7 tolerance such as 0.01 degrees. In the grid searching method for each (7, j-> ) anchor point 
needed per constraint, one fine parametric sweep of 7 is performed for a set constant . That sweep 
is then used to evaluate all of the constraints at that . 

A method to provide a good starting guess method for an iterative method, which has yet to be 
implemented, would be to use an iterative search method to locate a larger number of (7, jy) points. 
This method would start at a boundary point of the domain with an initial guess for a point on 
the curve, determined either by experience or from the undispersed database, and perform a 1-D 
line search to locate an initial point on the constraint. Then to find a subsequent point it would 
use a preset step size to traverse across the domain from one direction to the other based on the 
expected general orientation of the constraint in the domain (e.g. vertical, horizontal etc.). If the 
dependent data the line search is evaluating is continuous it may use gradient methods or even 
second gradient methods, while if it has discontinuities it will have to use zero-th order solvers 
such as the bisection method. After stepping in the expected direction of the constraint the routine 
would perform another 1-D line search to locate the (7, jy) on the constraint. Then the actual 
direction between the previous two points would be used as the step direction for the following 
step forward which would be followed by another 1-D line search. Storing the history of previous 
solutions should provide good initial guesses at each iteration step as most constraints can be closely 
modeled by low order polynomials. Essentially this is a method to walk along the constraint curve 
through the domain in order to characterize it. 


Computational Demands The computational demands are a function of the number of Monte 
Carlo anchor points, the number of Monte Carlo sets needed to identify all of the constraints, the 
number of constraints, the number of CPUs, and the required CPU time per Monte Carlo set. With 
respect to the grid search method, the iterative method may require a decreased number of Monte 
Carlo sets to find a single constraint, but may require an increase in the total number of of Monte 
Carlos to find all of the constraints when summed together. 

For the iterative method the total time would equate to that shown in Equation 9, assuming that the 
number of Monte Carlo iterations required to locate a constraint point is constant for each constraint. 


Ti me iterative method = MC Anchor Points • 


CPU time 
MC 


MC iterations 

• dim(g) ■ 

9i 


1 

No. CPUs 


(9) 


For the grid searching method the total time would equate to that shown in Equation 10. 

CPU time 1 

Time gr id method = MC Anchor Points • — — • MC sweep of 7 values 


MC 


No. CPUs 


( 10 ) 


Clearly from comparing Equation 9 and Equation 10 the condition where the run time crosses 
over would be that shown in Equation 1 1 . 


MC iterations 

9i 


■ dim(g) = Number of Monte Carlos in sweep of 7 values 


( 11 ) 


So if the iterative method converges on each of the constraints quickly, and the number of con- 
straints is small, it will reach a solution more quickly than the grid search method. The downfall of 
this method is that it requires a optimization routine, generates data that is not incremented spatially 
evenly, and does not provide insight into the sensitivity around each constraint limit in the domain 
space. 

Table 4 presents approaches to constraint curve fitting that offer combinations of point finding 
techniques, number of anchor points, and the order of the polynomial curve fit to apply. The beta 
entry in the first row of Table 4 represents using the already constructed undispersed database. In 
Table 4 four of the listed methods use the polynomial curve fit assessed in the undispersed database 
and anchors it to one Monte Carlo constraint point in ( 7 , ^) by translating it on the 7 axis, namely 
methods 1, 5, 6 , and 10. As such the constraint curve is really a pseudo-dispersed curve combining 
the speed and availability of the pre-existing undispersed database and the statistical accuracy of 
one Monte Carlo anchor point. 

For EM- 1 the entry corridor is often designed using the grid search approach labeled Method 7 in 
Table 4 with two high fidelity sweeps of 7 gathering enough information to characterize two points 
of each of the constraints. That is enough to fit a linear line. Suppose the needed number of Monte 
Carlo anchor points is two, the required number of runs per Monte Carlo is approximately 3400, 
the time per run is approximately 30 seconds, the number of CPUs available is approximately 90, 
and the number of Monte Carlos in 7 grid search is 200. The total time expectation for this process 
using Equation 10 is approximately 125 hours. On the other hand, cutting the needed number of 
Monte Carlo anchor points down to one, as suggested by Method 6 in Table 4 cuts the total time 
needed in half to roughly 62 hours. In that case the polynomial fit curve can be taken from the 
undispersed database and anchored (translated) by the determined anchor point. As yet, the idea of 
using more than two Monte Carlo anchor points, as in Method 8 and 9, has not been incorporated 
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Table 4. Constraint Curve Fitting Methods 


Method 

Point Finding 

Anchor Points 

Number of 

Polynomial Curve 

Polynomial Curve 

Technique 

Source 

Anchor Points 

Fit Source 

Fit Order 

Beta 

None 

Undispersed 

Database Specific 

Preset by User 

gi Specific 

1 

Iteratively 

Dispersed MCs 

1 

Undispersed Database 

gi Specific 

2 

Iteratively 

Dispersed MCs 

2 

Preset by User 

1 

3 

Iteratively 

Dispersed MCs 

3 

Preset by User 

2 

4 

Iteratively 

Dispersed MCs 

4 

Preset by User 

3 

5 

Iteratively 

Dispersed MCs 

Constraining 1 of n 

Undispersed Database 

gi Specific 

6 

Grid Search 

Dispersed MCs 

1 

Undispersed Database 

gi Specific 

7 

Grid Search 

Dispersed MCs 

2 

Preset by User 

1 

8 

Grid Search 

Dispersed MCs 

3 

Preset by User 

2 

9 

Grid Search 

Dispersed MCs 

4 

Preset by User 

3 

10 

Grid Search 

Dispersed MCs 

Constraining 1 of n 

Undispersed Database 

gi Specific 


due to their computational demands but that is planned for future work. For comparison consider if 
Method 1 was used which is the same as Method 6 except for using an iterative approach to finding 
anchor points. In that case, assuming there are 9 constraints and each constraint took consistently 
10 iterations apiece then the total time would be roughly 28 hours. 

The accuracy of the polynomial curve fits of the constraint curves, using the undispersed database, 
have been analyzed. The curves all appear to be modeled very well by polynomials of order 4 or 
3, effectively by order 2, and in the case of certain constraints even effectively by order 2 or 1. 
Polynomial approximations of constraints using order 5 and even some using order 4 suggested 
statistical over-fitting issues. The Table 5 shows a few example constraints with the deduced best 
order for a polynomial fit, along with the associated sum of Root-Mean-Square (RMS) to indicate 
the accuracy of the fits. 


Table 5. Constraint Curve Fitting: Guided, vi= 11.2, R=8890km (4800 n.mi) 


Constraint 

Side 

Polynomial Order 

Sum of RMS 

Guided Landing Accuracy 

Steep Side 

4 

0.0091057 

Aerothermal TPS Backshell Panel Temperature 

Top Side 

4 

0.0042080 

HSIR (e.g. Crew sensed acceleration) 

Steep Side 

3 

0.0064415 


As this mission is still early on in the design process, before the Preliminary Design Review 
(PDR), the trajectory design is still being evolved cyclically with feedback from other subsystems 
such as aerothermal, TPS and structures. The GN&C team provides the aerothermal team with 
sufficient trajectory data for it to assess the aerothermal flight conditions. Then it typically verifies 
that the aerothermal dependent variables being used to assess aerothermal requirements are accurate 
and sufficient from which to make vehicle level trades. Then if the designed Monte Carlo trajectory 
sets exceed the aerothermal teams limits the integration community must decide to either change 
the trajectory, change the aerothermal modeling approach, or to change the vehicle hardware such 
as the TPS, gap fillers, outer mold line. 

After the previous cycle the aerothermal team determined the following set to determine TPS 
sizing drivers. In previous design cycles GN&C has delivered trajectory data to the aerothermal 
team to assess for a large Monte Carlo run matrix composed of TPS size driving combinations of 
the following: or 'y EI:Steeplirnit , guidance mode, El inertial velocity magnitude, El 
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range, jj, j EI , and landing site. That cycle identified a sensitivity to j EI that is driving the TPS 
sizing with decreasing r ) EI (steepening) leading to increasing TPS mass, and conversely increasing 
7 b/ (shallowing) leading to decreasing TPS mass. 


ENTRY CORRIDORS 


After a recent cycle of creating and evaluating integrated entry corridors certain system level as- 
sumptions were made to permit decrease ^ values and increase 7 m i n to decrease mass required 

for CG ballast placement and TPS sizing. In Figure 3 two pairs of entry corridors 73 Mach =25 vs - 
7 b/ are shown for two different entry ranges that are 4630 km and 8890 km. 



Inertial Topocentnc Flight Path at El (d) 



(a) Entry Corridor for Guided, v/=10.95 ^2, R=4630 km (b) Entry Corridor for Guided, u/=10.95 R=4630 km 

(2500 n.mi.). Neglecting SM Debris Disposal Constraint (2500 n.mi.) 



Inertial Topocentnc Flight Path at El (d) 
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Inertial Topocentnc Right Path at El (d) 


1 

< 

158 6 o 
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(c) Entry Corridor for Guided, tn=10.95 ^ , R=8890 km (d) Entry Corridor for Guided, u/=10.95 R=8890 km 

(4800 n.mi.). Neglecting SM Debris Disposal Constraint (4800 n.mi.) 


Figure 3. Entry Corridors 


In the plots in Figure 3 the constraints have been represented by areas of color compressing the 
flyable corridor colored white. The constraints are identified in the legend along with the calculated 
corridor widths at three values. After doing more in-depth analyses of the low values in Figure 
3(a) and Figure 3(b) the decision was made to hold the minimum at 0.23 as indicated in Table 3. 
The program has decided to neglect the SM debris disposal limit as it impacts the corridor design. 
To assess the impacts of this choice see Figure 3. The SM debris disposal line is indicated by the 
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blue line. The margin that the program has chosen to cover for r ) EI corridor width of 0.15 degrees 
is indicated as an offset red dashed line from the purple line representing the Ballistic 2er lofting 
constraint. Hence, when the red dashed line (the 0.15 degree margin) is steeper than the SM debris 
disposal constraint the margin is adequate, but when it is shallower than the SM debris disposal 
constraint line the margin is inadequate and the design needs to move to the real SM boundary. The 
Monte Carlo anchor points that were used to anchor the constraint polynomials from the undispersed 
database are shown as squares. The squares that are blackened indicate the four extreme cases in the 
corridor that GN&C assessed and delivered to other subsystems for further analysis. It is clear from 
each pair of figures, Figure 3(a) with Figure 3(b) and Figure 3(c) with Figure 3(d), above that the 
assumption to ignore the impacts of SM debris disposal is a serious risk. To deal with this risk EM- 1 
is planning to perform a CM raise maneuver that it was not planning to previously. The combined 
CM and SM vehicle will target a steeper than nominal 7 S/ , separate from each other leaving the 
SM at that return 'y EI , and then the CM will use its RCS jets to burn to shallow up its own r y EI . This 
maneuver will ensure that the SM is not a risk to the public and that the CM can target a shallower 
than it otherwise would be able. 


CONCLUSION 

The characterization of a spacecraft’s entry corridor is governed by coupled nonlinear differential 
equations, active control, and inputs from many other subsystems. The characterization is sensitive 
to changes from other sub-systems such as mass properties, propulsion, aerodynamics, TPS, struc- 
tures, etc.. Depending on the magnitude and direction of the changes the associated spacecraft entry 
corridor will change in different ways. Reacting repetitively to new inputs has lead to a large amount 
of strategic automation leveraging pre-existing understanding to characterize corridors quicker and 
more definitively each iteration. This improvement has lead to the entry corridor being the one of 
the first steps in assessing a spacecraft’s re-entry capabilities. 


FUTURE WORK 

The authors acknowledge this work is in support of a mission that has yet to reach Preliminary 
Design Review (PDR) and will surely make design changes in the coming design cycles that cannot 
be foreseen and taken into account. Thus this work will be continued in order to improve the robust 
capability and speed of evaluating re-entry flight corridors for vehicle design changes. It is the 
intent that as the design process advances other trajectory requirements will be taken in account in 
the entry corridor design including drogue parachute deployment flight condition, main parachute 
deployment flight condition, landing and touchdown heading angle among others. 
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NOTATION 


c Flight Path Angle Corridor Width 

Lift-over- Drag Ratio at Mach Number of 25 of Vehicle 
m Vehicle Mass at Entry Interface 
R Downrange 

v Earth Centered Inertial Velocity Vector Magnitude of Vehicle 

6 Geodetic Latitude 

7 Inertial Topocentric Flight Path Angle 
-0 Inertial Topocentric Azimuth 

A Longitude 
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